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Anomalous diffusion mediated by atom deposition into a porous substrate. 
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Constant flux atom deposition into a porous medium is shown to generate a dense overlayer and 
a diffusion profile. Scaling analysis shows that the overlayer acts as a dynamic control for atomic 
diffusion in the porous substrate. This is modeled by generalizing the porous diffusion equation 
with a time-dependent diffusion coefficient equivalent to a nonlinear rescaling of time. 
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Thin film depositions on substrates are important in 
many physical processes and applications. Moreover, de- 
position on porous substrates is particularly useful for 
catalytic systems [H in which atomic deposition carried 
out by plasma sputtering is then coupled with the diffu- 
sion of the atoms into the porous substrate. Diffusion 
in a porous medium is generally anomalous and is char- 
acterized by the mean square displacement < z 2 (t) >: 
it evolves with a power law in time t a different from 
the well known linear behavior ( a = 1) for normal dif- 
fusion. Anomalous diffusion [1, 0| is a general process 
that can be observed in many domains such as transport 
in porous [|[ and/or fractal media [f|, surface growth 
0, solid surface diffusion § or hydrodynamics (rotat- 
ing flows Q, turbulence [lol ] or diffusion in an array of 
convection rolls [111 H^). It can be described by differ- 
ent models that involve space-dependent diffusion coeffi- 
cients [l,ll,[ll,[l3,[H,[H,Il3. In these models anomalous 
diffusion is included through the propagator (the solution 
of the diffusion equation starting with a Dirac distribu- 
tion at t — 0) that exhibits both power laws and stretched 
exponential behaviors. Here, we study platinum deposi- 
tion by plasma sputtering on a porous carbon substrate. 
The platinum atoms are sputtered by the plasma ions 
and travel to the substrate so that both deposition and 
diffusion processes can depend on plasma operating pa- 
rameters. The experimental results consist in a measure 
of the density profiles of the deposited matter at differ- 
ent times. The goal of this Letter is to characterize this 
diffusion-deposition process in a porous carbon medium 
using the time evolution of the platinum density profile. 

Platinum atoms are deposited by plasma sputtering 
into a porous carbon layer supported (See Fig. [I} on 
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a carbon cloth. This porous layer is a few tens of mi- 
crons thick and is composed of randomly stacked car- 
bon nanoparticles (Vulcan XC 72) and PTFE particles 
brushed onto the carbon cloth. The specific area is 
15 m 2 g _1 before Pt deposition and slightly lower at 13 
m 2 g _1 after Pt deposition. Examination of FigfT] shows 
the pore size reduction. Platinum atoms are deposited 
onto this porous carbon layer and diffuse in the course 
of deposition. Theplatinum deposition reactor is a pre- 
viously described [If, plasma sputter system which de- 
livers a constant deposition rate. An argon plasma is 
created in the stainless steel deposition chamber 18 cm 
inner-diameter and 25 cm-long by using an external pla- 
nar antenna (also known as TCP antenna). The porous 
layer is placed on a movable grounded substrate holder 
in front of the platinum sputter target with a target- 
substrate distance of 4.5 cm. The porous layer is thus 
exposed to a flux of sputtered platinum, into which it 
diffuses under the operating conditions. The diffusion 
process was studied at argon pressures of 0.5 and 5 Pa, 
leading to a mean kinetic energy Ek of the sputtered 
Pt atom around 7 eV Q and 0.04 eV respectively [ls| . 
The resulting platinum depth profiles were measured us- 
ing Rutherford Backscattering Spectroscopy, which gives 
an indirect measure of the average density as a func- 
tion of the depth z. The density is obtained by fitting 
the experimental spectrum with a spectrum derived from 
a defined profile function. Firstly, we observe that the 
platinum profile can in fact be decomposed in two dis- 
tinct regions: a growing layer above the porous medium, 
of mean thickness h(t) — zo(t) and constant density 
p(z,t) = Zi(t),—zo(t) < z < and a density profile 
p(z, t), z > in the porous medium (which extends in the 
z > domain, z = corresponding to the porous inter- 
face) . The Pt depth profile is thus deduced by minimizing 
the difference between the experimental RBS spectrum 
and the simulated RBS spectrum using the profile de- 
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FIG. 1: Scanning Electron Microscopy top view of porous 
carbon layer a) before Pt deposition b) after Pt deposition 
(Courtesy D. Cot, IEMM CNRS-ENSCM-Univ. Montpellier 
W 



fined by Eq. [TJ It corresponds to the known solution of 
anomalous diffusion processes @ i.e. a stretched expo- 
nential. 



p{z,t) = Zi(t), 

2 + 

p(z,t) = Zi(t)e~*5W, 



-Zq < z < 
z > 



(1) 



where 9 is the dimensionless coefficient which character- 
izes the anomalous diffusion behavior [ID, [lj, fl5| .This 
fitting procedure leads to the error bar of the exponent 9 
as reported hereafter. The density profile is continuous at 
the interface between the growing platinum layer and the 
porous medium at z = so that we have p(0,t) — Z\(t). 
Z^ify is the time dependent spreading of the stretched ex- 
ponential. Such a profile suggests that the particle flux 
delivered by the plasma sputtering cannot be absorbed 
directly by the porous layer with the result that a frac- 
tion of the incident atoms have to deposit at the interface. 
It also implies that the diffusion process into the porous 
layer is a consequence of both the diffusion of the plat- 
inum atoms that have penetrated into the substrate and 
the diffusion of those absorbed at the interface. 
Some fitted platinum depth profiles at different deposi- 
tion times t in the porous medium are displayed in Fig. 
H The best fits correspond to 9 = -1.45 ± 0.05 « -3/2 
which suggests a super-diffusive behavior in the porous 
medium (see below). Moreover, given 9 — —1.5, Z\(t) 
and Z2{t) asymptotically follow power law behaviors with 
time (Z x ~ t m and Z 2 ~ i p ) with m being 0.40 ± 0.05 
and p being 0.20 ± 0.05 as shown in Fig. EJa). h(t) 
also obeys a less well defined power law behavior with 




-100 



100 200 300 

depth z (nm) 



400 



500 



FIG. 2: (Color online) Pt depth profiles in porous carbon at 
different successive deposition times. z = corresponds to 
the porous carbon layer surface 



~ t n ,n = 1.1 ± 0.1 . The mass of platinum inside the 
porous medium is found to scale as i - 8 * - 1 while the Pt 
mass of the surface overlayer, Z\(t).h(t), scales as t 1 - 5 ^ - 1 
as shown on Fig. [3Jb). However, we obtain an almost 
linear evolution of the total mass (~ t 11±01 ) of platinum 
deposited in and on the porous medium, as expected by 
the experimental process. At the higher pressure, 5 Pa, 
9 = -1.33 ± 0.05,ra = p = 0.25 ± 0.05 and n = 1.0 ± 0.1. 
Hence, a change of deposition conditions yields a change 
of Pt concentration profile within the porous carbon. To 
help explain our results, we use a classical model of the 
anomalous diffusion process [HI, 0, EH, EH in fractal 
and porous media, considered here in one dimensional 
space accounting for the in plane averaging (we can dis- 
card nonlinear diffusion in the porous media from the 
experimental results) resulting in the following general- 
ized diffusion equation for the density profile p(z 1 t) for 
z > 0: 



dp 

at 



d_ 

dz 



Kodp 

z e dz 



(2) 



where 9 is a real parameter characterizing anomalous dif- 
fusion and is a priori unknown. This equation is obtained 
by postulating a distance dependence of the diffusion 
coefficient in Fick's law. Classical Brownian diffusion 
holds for 9 — 0, whereas positive (negative) 9 corresponds 
to sub(super)-diffusive dynamics where the mean-square 
displacement scales as < z 2 >oc t 2 /( 2+e ). The general 
solution for this equation (the so-called propagator) ex- 
hibits a power law and stretched exponential behavior 
and is given by @, O, M, M, [H : 



p(z,t) oc K (2 + 9) t 



-1/(2+6) 



exp 



.,2+0 



K (2 + 9)H 



(3) 

which holds for a normalization condition consistent with 
mass conservation. Since the deposition is carried out 
with a constant flux, the total mass is expected a priori 
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FIG. 3: (Color online) For the 0.5 Pa deposition, a) thickness 
of the Pt layer h(t) above the porous medium and value of 
Z\(t) and Ziif) from the fitted law (fJJ) as a function of time 
(log- log plot). The lines correspond to the power laws Z\ (t) ~ 
t 0A ,Z 2 (t) ~ t 2 and h(t) ~ t 11 ; b) Amount of Pt mass in the 
porous medium calculated using (fj} as a function of time (log- 
log plot). 



to grow linearly and the solution to our problem should 
obey the boundary condition at z = 0: 



^l2e = Co 

z dz 



(4) 



where Co is the constant flux of platinum imposed by 
the plasma sputtering. Therefore, the solution can be 
obtained using the propagator solution ([3J through a lin- 
ear superposition consistent with the boundary condition 
([3]). It is in fact more convenient to investigate the solu- 
tion in the following self-similar form: 



pM) = f7(£) 



(5) 



While the propagator obeys a + (3 = (constant mass so- 
lution), the constant flux solution must have a + /3 = 1. 
Inserting ([5]) into the diffusion equation we obtain 
a = while (3 is fixed by the flux condition. The 
self-similar function /(£) is then the solution of the ordi- 
nary differential equation (£ = z/t a being the self-similar 



Pf-<xU' = ^{U"-6f) 



(6) 



which gives the propagator solution ([3]) for a + (3 = 0. 
Moreover, it can be seen that the density profile for 
large z follows the stretched exponential law ([I]) with 
Zi(t) tx and Z 2 (t) = K (2 + 9) 2 t. 
The predictions of this simple model do not agree with 
the experimental results: from 9 = —1.5, one would ob- 
tain a = 2 and (3 = — 1 since we have the condition 
a + /3 = 1. Hence, there is a marked discrepancy between 
the experimental measurements (a ~ 0.4 and (3 ~ 0.4) 
and the model even when the error bars of the exper- 
iments are taken into account. However, the underly- 
ing physical processes of atomic diffusion and the mea- 
sured density profiles suggest that the general diffusion 
equation proposed here is a good framework for model- 
ing atomic deposition through plasma sputtering. It is 
thus tempting to model the effect of the growing external 
platinum layer and of the plasma by keeping the same 
equation with the anomalous diffusion coefficient K(t) 
and a flux of mass at z = 0, both time dependent. Al- 
though we are far from a detailed microscopic derivation 
of the model, it is postulated here that the growth of the 
platinum layer has a screening effect on the diffusion and 
thus alters the mass flux towards the porous medium. In 
particular the platinum overlayer acts as a reservoir for 
diffusion. In addition, the temperature distribution, not 
accounted for directly by this model, is certainly influ- 
enced by the presence of the platinum layer. Since the 
diffusion coefficient is a function of the temperature, it is 
probably time dependent. This time dependent diffusion 
coefficient can also be understood as a time rescaling and 
it is needed in the framework of the diffusive equation 
since it is the only way to change the time dependence 
in the exponential law ([3]). We thus propose to model 
the deposition/diffusion dynamics through the previous 
equation Q with a time dependent diffusion coefficient 
K{t): 



dp 

at 



d_ /K(t) dp 
dz \ z e dz 



(7) 



The presence of the external growing layer is also ac- 
counted for by an a priori time dependent boundary con- 
dition: 



K(t) dp 
z 6 dz 



C{t) 



(8) 



To illustrate the results with no loss of generality, we 
will seek a power law behavior for a sufficiently time. We 
introduce two additional exponents e and 7 such that the 
flux of atoms and the diffusion coefficient follow: 

C(t) = C t £ K{t) = Kof 1 (9) 

The relations between the exponents are straightforward: 

7+1 



and a + f3 — 1 + e 



(10) 
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Eq. ([7]) can be understood by rescaling the time following 
t = t 1+1 in terms of the initial diffusion equation @: 



dp 
8t 



d ( K o 1 dp 



dz V 1 + 7 z e dz 



(11) 



And we obtain the stretched exponential behavior for 
large z in agreement with the experimental measure- 
ments: 



p(z, t) oc 1r exp 



(1+7) 



K (2 + 9)H x +i 



(12) 



Within this general approach, we can now describe our 
experimental results. For simplicity's sake and to illus- 
trate the model with no loss of generality we estimate 



from the experimental fits that 9 = 
j3 = 0.4 which leads to: 



= 0.4 and 



7 = -0.8 



and 



e = -0.2 



This corresponds to a time decreasing diffusion coefficient 
such as t~ - 8 and of the flux at z = such as t~°- 2 . The 
diffusion coefficient is expected to decrease versus time, 
while the porous medium is gradually filled by the plat- 
inum atoms. This leads to a more difficult diffusion. At 
the z = interface, the growing external layer provides 
an additional source for the diffusion of atoms inside the 
porous substrate. Moreover a decreasing flux indicates 
that the growing external layer is gradually screening the 
influence of the plasma, and that diffusion becomes less 
efficient due to pore filling as mentioned above. These 



effects arc also observed at 5 Pa, with 9 = 



a = 0.5 



and (3 = 0.25 which leads to: 7 = -0.83 and e = -0.25. 
These last two scaling exponents are very close to those 



at 0.5 Pa. Consequently, the self similar exponents de- 
pend on the deposition conditions whereas the critical 
exponents of the diffusion are unchanged. This means 
that the diffusion equations (f7ll8 TTTj) are robust enough 
to provide a general frame for deposition/diffusion mech- 
anisms. 

On the basis of experimental results and a general dif- 
fusion model, we have proposed a new scenario for mass 
deposition and diffusion on a porous substrate. As the 
flux of sputtered atoms cannot be absorbed immediately 
by the substrate, an overlayer grows. The combined ef- 
fects of the deposition that change the porosity properties 
of the substrate and of this growing layer that can act as 
a reservoir for diffusion can be accounted for by general- 
izing a standard model of diffusion in a porous medium 
with constant diffusive properties (9 is constant). This 
allows time dependent coefficients to be introduced. This 
model provides a clear understanding of the experimental 
results: in particular, we show that the diffusion coeffi- 
cient in the porous substrate and the flux at the inter- 
face decrease with time due to the atom deposition in 
the substrate. This general model can be also be applied 
to many other diffusion processes where the inner porous 
structure of a material is modified by a penetrating flux 
of matter, as for example, for filtration through porous 
membranes, sedimentation processes, ... 
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